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Abstract 

We investigate the dependence of the time delays for the large-separation gravi- 
tationally lensed quasar SDSS J1004+4112 on the inner mass profile of the lensing 
cluster. Adopting the mass model whose innermost density profile is parameterized as 



O 



X 



p oc r a , we derive a series of mass models which can fit observational data and then 



compute the probability distribution functions of time delays. We find that larger a 
q has longer time delays, longer tails at the higher end of the probability distribution, 

■ and larger model uncertainties. The ratios of time delays slightly depend on the slope 

a. Among others, time delays between images C and A (or B) have little dependence 
on the inner slope, particularly when the time delays are short. The dependence of 
time delays on a is well fitted by a linear form, which reflects well-known degeneracy 
between the mass profile and time delays. We perform a Monte-Carlo simulation to 
illustrate how well the inner slope can be constrained from measurements of time 
delays. We find that measurements of more than one time delays result in reasonably 
tight constraints on the inner slope (a a ^ 0.25), while only one time delay cannot 
determine the inner slope very well. Our result indicates that time delays indeed 
serve as a powerful tool to determine the mass profile, despite the complexity of the 
lensing cluster. 

Key words: cosmology: theory — galaxies: halos — galaxies: clusters: general 
- dark matter — galaxies: quasars: individual (SDSS J1004+4112) — gravitational 
lensing 
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1. Introduction 

Recent high-resolution iV-body simulations in the Cold Dark Matter (CDM) universe 
have suggested that dark matter halos are described by a universal mass profile (Navarro et 
al. 1996, 1997, hereafter NFW). Higher-resolution simulations have revealed that the inner- 
most region may have steeper inner slopes than originally suggested and may not be uni- 
versal (Fukushige, Makino 1997, 2001, 2003; Moore et al. 1999; Ghigna et al. 2000; Jing, 
Suto 2000; Power et al. 2003; Fukushige et al. 2004; Hayashi et al. 2004; Diemand et al. 
2004; Tasitsiomi et al. 2004; Reed et al. 2005). While the accurate value of the inner slope 
and the universality remain controversial, the existence of deep potential well at the center of 
dark matter halos in the CDM model appears quite robust. Here the important fact is that the 
inner profile of dark halos is extremely sensitive to nature of dark matter: For instance, if one 
introduces self-interacting cross sections of dark matter particles, the central density profile 
becomes much shallower than expected in the standard CDM universe (Spergel, Steinhardt 
2000). Therefore it serves as a powerful test of collisionless CDM paradigm. Such observa- 
tional studies to constrain parameters of the density profile have been conducted. For instance, 
assuming the hydrostatic equilibrium the mass profile of a dark matter halo can be evaluated 
from the gas density and the temperature profiles of intra-cluster medium (e.g., Sato et al. 
2000); weak lensing reconstruction can reproduce the mass profile of a cluster from ellipticities 
of galaxies behind the cluster (Kaiser, Squire 1993; Schneider 2005). In particular, rotation 
curve observations of dark-matter-dominated dwarf and low surface brightness disk galaxies 
have indicated that they favor mass profiles with a flat density core, being inconsistent with 
the NFW profile (Salucci, Burkert 2000; de Blok, Bosma 2002; Swaters et al. 2003; Gentile et 
al. 2004). However, this apparent discrepancy might be because of bias in gas rotation speed 
(Hayashi et al. 2004) or the destruction of central cusps by the formation of a primordial bar 
(Weinberg, Katz 2002). Thus, investigations of the halo density profiles in different systems 
has been one of the most essential issues. 

Among others, one of the most promising methods to test the NFW profile would be 
strong lensing by clusters of galaxies. Gravitational lensing is unique in the sense that it allows 
one to detect the mass distribution directly (Refsdal 1964; Bourassa et al. 1973, 1975, 1976; 
Wambsganss, Paczyhski 1994; Keeton et al. 1997; Mao, Schneider 1998; Chiba 2002; Kawano 
et al. 2004; Kochanek, Dalai 2004). In particular, strong lensing by clusters is a indispensable 
tool to probe the innermost region of dark matter halos, since most baryons in clusters remain 
hot and diffuse and therefore the density profile of clusters can be well approximated by that 
of dark halos seen in iV-body simulations. For instance, giant arcs of background galaxies 
with different redshifts (Einstein radii) can strongly constrain the mass profiles of clusters 
(Broadhurst et al. 2005). This is because the mass inside the Einstein ring with radius 
at which arcs build up is proportional to r| and thus multiple arcs with different redshifts 
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determine the masses within the different radii, equivalently the mass slope. Multiply imaged 
quasars due to clusters of galaxies offer another direct probe of the mass profiles of lensing 
clusters. Indeed, they have several advantages over giant arcs. First, for lensed quasar systems it 
is easier to correct selection bias of lensing clusters (i.e., difference between lensed and unlensed 
cluster populations; see e.g., Hennawi et al. 2005) because of well-known source population of 
quasars. Second, we can measure time delays for lensed quasars, which serves as additional 
strong constraints on the mass distribution. 

Recently, the first example of such quasar-cluster lens system, SDSS J1004+4112, has 
been discovered (Inada et al. 2003; Oguri et al. 2004) in the Sloan Digital Sky Survey (SDSS; 
York et al. 2000). It has an unusual separation of ~ 15 arcsec, which is more than twice larger 
than the second largest lens Q 0957+561 (Walsh et al. 1979). The quasar and lensing cluster 
have redshifts of 1.74 and 0.68, respectively. The central region of the lens cluster is dominated 
by the brightest cluster galaxy. The fifth image of the lensed quasar, which constrains inner 
mass distribution of the the brightest cluster galaxy strongly, was discovered by Inada et al. 
(2005). 

How did the lens system SDSS J1004+4112 constrain the inner mass distribution of 
the dark matter halo? Oguri et al. (2004) modeled the lens with singular isothermal ellipsoid 
(SIED) plus NFW, and found that the image configurations and fluxes are reproduced well. 
More generally, Williams, Saha (2004) adopted a free-form mass reconstruction technique (Saha, 
Williams 1997) and pointed out that the quadruple images alone do not give useful constraints 
on the inner slope of the dark matter halo. Therefore, we need additional information: One 
such information comes from arcs of background galaxies which have already observed by 
Sharon et al. (2005). Arcs at several radii (redshifts) can constrain the inner slope since we 
can estimate the mass inside a arc from the the arc position and shape. Another information, 
which is unique for lensed quasar systems, is time delays between each images. It has been 
shown that the Hubble constant and the slope of the density profile are both sensitive to 
time delays (Wambsganss, Paczyhski 1994; Oguri et al. 2002; Wucknitz 2002; Kochanek 2002; 
Oguri, Kawano 2003; Kochanek et al. 2005), therefore by assuming the value of the global 
Hubble constant, which is now determined with better than 10% accuracy (Freedman et al. 
2001; Spergel et al. 2003), we will be able to obtain useful information on the mass profile. 
The image configuration of the system is very similar to that of PG 1115+080 (Weymann et 
al. 1980), however the separations are scaled by a factor of ~ 8 and hence the time delays are 
longer by a factor of ~ 8 2 , which means that the shortest and longest time delays would be ~ 10 
days and ~ 1000 days, respectively. However, the complexity of the cluster mass distribution 
makes the possible range of time delays quite large (Oguri et al. 2004; Williams, Saha 2004). 
Therefore it is not obvious whether time delays can really give meaningful constraints on the 
inner profile of the dark matter halo. In this paper, we investigate in detail the capability of 
measuring time delays and constraining the inner mass slope of the cluster. 
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This paper is organized as follows. In section 2 we describe our model to predict time 
delays for SDSS J1004+4112. Section 3 gives the predictions of time delays and the probability 
distributions of them. Then, a likelihood analysis is performed in section 4 to see how time 
delays constrain the radial slope of the cluster mass. We finally discuss the results and give 
conclusions in section 5. 

2. Lens Mass Modeling 

A source at y and an z-th image at Xj are related through the lens equation 

y = x,-V^(x i ), (1) 

where ipfc) is the projected lens potential (Schneider et al. 1992). The vectors are defined on 
the sky (lens and source planes). The lens potential can be expressed as 

V>(x) = ~ /dV«(x')ln|x-x'|, (2) 

7T J 

where re(xj) is the dimensionless surface mass density, so called convergence. More specifically, 
the convergence is proportional to the surface mass density S(x) as 

= -ir 1 , (3) 

^cr 

where S cr is the critical surface mass density defined by 
c 2 D 

Scr = 1. (4) 

with D d , D 8 and D ds being the angular diameter distances to the lens, to the source, and from 
the lens to the source, respectively. The gravitational lenses magnify the flux of i-th image by 
the factor of /i(x;): 

-i 



/i(Xj 




(5) 



(6) 



The relative time delay between z-th and j-th images is then calculated from 

A 1 + z A DgDs 
c D ds 

where Axj = x, — y. Provided that one knows redshifts of the source and the lens, position 
and flux of the source, and mass distribution (potential) of the lens as well as the cosmological 
parameters, one can calculate observable quantities such as the number of images, image posi- 
tions, fluxes, and time delays. Turning the problem around, by measuring the redshifts, image 
positions and fluxes, and relative time delays, one can constrain the mass distribution of the 
lens. 

In this paper, we consider a specific quadruple lens system, SDSS J1004+4112, the first 
quasar multiply lensed by a central part of a massive cluster. There has been a variety of 
studies to investigate the structures of the lensing cluster. Oguri et al. (2004) firstly performed 
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an enormous amount of modeling of the system with realistic two-component models. Provided 
that the galaxy has a SIED and the cluster has an elliptical version of NFW profile with an 
external shear, respectively, they concluded that: 1) there is a offset between the centers 
of the brightest galaxy and the cluster; 2) a wide range of the models can reproduce the 
position angle of the galaxy well; 3) the elongation of modeled cluster are also similar to 
the observed distribution of cluster galaxies; 4) they found a large tidal shear (~ 0.2) which 
suggests significant substructure in the cluster; 5) there is an enormous uncertainty in the 
predicted time delays between the quasar images, which indicates that measuring the delays 
would greatly improve constraints on the models; 6) measuring the time delay between image 
A and B would determine the temporal ordering such as C-B-A-D and give an expected value 
of the delay between image C and D; 7) the scale length of the dark matter of the cluster would 
be r s > 30" from predicted relation between r s and the lensing strength (see below for the 
definition of r s ). On the other hand, Williams, Saha (2004) presented free- form reconstructions 
of the lens with constraints of the image positions and physical conditions. The modeling gave 
some important results: 1) the projected cluster mass profile is consistent with being r - - 3 - -°- 5 ; 
which can be fitted with either the NFW or a flat core model; 2) the residual mass maps created 
by subtracting the circularly averaged surface mass density shows the significant substructures; 
3) D-A-B-C time ordering results in shift ~ 3" of the lens center and then tends to allow a lot of 
models with spurious extra images — the time ordering seems very unlikely; 4) a measurement 
of time delay between A and B would serve as a test of the shallow mass profile. Inada et al. 
(2005) discovered the fifth central image and tested the ability of the image to constrain the 
mass profile of the bright cluster galaxy. They assumed a power-law density profile p(r) oc r~ 7 
which potential was set to have an elliptical symmetry for the galaxy and concluded that 
the central bright cluster galaxy cannot have steeper mass profile than isothermal (7 < 2) to 
reproduce the flux. Sharon et al. (2005) discovered multiply imaged arcs of three galaxies at 
high redshifts. Their different redshifts correspond to different radii of Einstein rings because 
of different critical surface mass densities. Since multiply imaged arcs constrain the mass inside 
the arcs, they obtained the average projected surface mass slope consistent with NFW (~ —0.5) 
from their simple (preliminary) one-component modeling of the dark matter halo. We note that 
the lensing probability can constrain the mass profile as well (Oguri et al. 2004; Oguri, Keeton 
2004), but the constraint degenerates strongly with cosmological parameters. 

Following Oguri et al. (2004), in this paper we consider the two-component model that 
consists of the brightest cluster galaxy and dark matter components. However, an important 
extension of the model in this paper is that we allow various inner slopes of dark matter 
components to study their effects on time delays. More specifically, we model the brightest 
cluster galaxy, which Oguri et al. (2004) called Gl, as an SIED 

K (%,y) = W7, (7) 
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Table 1. Constraints on Mass Models 



Object 


x (arcsec)* 


y (arcsec) * 


Flux (arbitrary) ^ P.A. (deg) * 


A 


0.000 ±0.001 


0.000 ±0.001 


1.000 ±0.200 


B 


-1.317±0.002 


3.532 ±0.002 


0.732 ±0.146 


C 


11.039 ±0.002 


-4.492 ±0.002 


0.346 ±0.069 


D 


8.399 ±0.004 


9.707 ±0.004 


0.207±0.041 


Gl 


7.114±0.030 


4.409 ±0.030 


-19.9 ±20.0 



* TEe positive directions ol x and y arc denned by West and North, respectively, 
t Errors are broadened to 20% to account for possible systematic effects, 
t Degree measured east of north. 



where £ = yx 2 + y 2 /q 2 , b is a characteristic deflection angle which is related to the velocity 
dispersion a v of the galaxy by 

b = ^M, (8) 



A V c 

and q is an axis ratio of the projected mass distribution. The separation of the lensed images 
is so large that the small deviation from the isothermal profile have little effect on predictions 
of the time delays. The cluster of galaxies is modeled as a generalized NFW profile (GNFW) 

n = Ps (g) 

where r s is a scale length, p s is a characteristic density (computed from the scale length r s and 
the virial mass M vir ), and a is an inner density slope (p oc r~ a ). The surface mass density of this 
GNFW profile is obtained by integrating equation (9) over z. As done by Oguri et al. (2004), 
an external shear is added in order to approximately include the effects of possible complex 
structure in the outer region of the cluster. 

In order to exclude unphysical situations, we assume that b is smaller than 2.25 arcsec 
because any galaxies should have the velocity dispersion smaller than 400kms" 1 . The inner mass 
slope a is restricted to 0.5, 0.75, 1.0, 1.25, 1.5 mainly because of the computational limitation. 
We note that NFW originally proposed a = 1 universal model and recent N-body simulations 
have suggested that a is approximately between 1 and 1.5. We fix r s = 40", because Oguri et al. 
(2004) concluded that r s > 30", which is much larger than the values of the distances between 
the lens center and the images. In an inner region (r <Cr s ), profiles with different scale lengths 
are not much different from each other while varying a has much larger effects. This is clearly 
demonstrated in Figure 1, where we show the predicted time delays of r s = 40", 60" (each point 
in the Figure corresponds to one of statistically acceptable models; see below for more details). 
This figure implies that the dependence on the scale length is not important in the case of SDSS 
J1004±4112, as long as the scale length is reasonably large, r s > 30". Therefore, the number of 
the model parameters for each a is 15: the position, mass, ellipticity, and position angle of the 
central galaxy Gl; the position, mass, ellipticity, and position angle of the cluster (GNFW); 
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Fig. 1. Predictions of the three time delays for models with different scale length. We fix the inner density 
slope to a = 1. Green crosses show predictions for r s = 40" and black squares denote those for r s = 60". 

the amplitude and position angle of the external shear; the position and flux of the source. 

For observed positions and fluxes of the quasar and the central galaxy, we use the data 
of Inada et al. (2005), which is summarized in Table 1. For the purpose of investigating relative 
time delays between the large-separation images, information of the fifth image is not included 
in our calculation, because its position and flux are supposed to be governed mainly by the 
detailed mass distribution of the bright cluster galaxy Gl rather than that of the dark matter 
halo. Actually, we assume the galaxy to have an isothermal mass profile and then our profile 
has passed the test of Inada et al. (2005). In addition to image positions and fluxes of the the 
quasar and the galaxy, the data of a position angle of the galaxy (9 g ~ —19°. 9 ± 20°. 0) is used 
to constrain our model since the projected mass and light are generally aligned to each other 
(Keeton et al. 1998). As a result, there are 15 constraints on mass models. 

Since the numbers of the model parameters and the observables are equal to each other, it 
is expected that a wide range of parameter space can fit the data (e.g., Keeton, Winn 2003; Oguri 
et al. 2004). Following their work, we derive a series of acceptable models as follows. First, we 
start from a random starting point in the parameter space. We put the parameters initially 
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Fig. 2. The A% 2 distribution as a function of the inner slope a. The minimum % 2 was achieved at a = 1.5, 
and the value was x 2 = 0.01. 

on the range 0" < b < 2". 25, 6". 964 < x(Gl) < 7".264, 4".259 < y{Gl) < 4".559, < e g < 0.9, 
-59°. 9 <0 g < 20°. 1, < k s = p s r s /S cr < 1.5, < e c < 0.9, -90°. < 6 C < 90°. 0, < 7 < 0.8, 
— 90°. < # 7 < 90°. 0, where e g is the ellipticity of the galaxy, e c and 9 C are the ellipticity and 
its position angle of the cluster, and 7 and # 7 are the amplitude and its position angle of the 
external shear, respectively. We then perform a x 2 minimization and find a local minimum in 
the x 2 surface. We adopt lensmodel package (Keeton 2001) to solve the lens equation and to 
perform x 2 minimization. If the x 2 °f the minimum is < 11.8, we regard it as an acceptable 
model, and pick it up 1 . Any model with b > 2". 25 or unusually large ellipticities e = 1 — q > 0.9 
is excluded so as not to include unphysical models in our analysis. Because of the result 3) of 
Williams, Saha (2004) described above, we allow only the models that predict C-B-A-D time 
ordering. By repeating this process from random numerous starting points, we obtain a group of 
acceptable models. Most (> 90%) of local minima in the parameter space are excluded because 
they are not physical or give bad fits statistically. For each a, we determine 100 acceptable 
models accordingly. 

Throughout the paper we adopt a flat lambda-dominated universe with (Q Q , Q\) = 
(0.3,0.7), where Qq is the density parameter of matter and Qa is the dimensionless cosmo- 
logical constant. However, our predicted time delays can be converted to those in any other 
cosmologies via equation (6). The Hubble constant in units of 70 km s _1 Mpc -1 is denoted by 
h 70 . 

1 The value of 11.8 represents the 3<r limit in a projected two-dimentional parameter space, when Xmin = 0- 
We note that the value was also adopted in Oguri et al. 2004. 
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Fig. 3. The A% 2 distributions of our 100 acceptable models for each a. 
3. Predictions of Time Delays 

Using the method described in the previous section, we fit the gravitational lens system 
SDSS J1004+4112 with varying inner mass slopes, a = 0.5, 0.75, 1.0, 1.25, 1.5. First, we regard 
a as a parameter and compute A% 2 = \ 2 ~ Xmm- For each a, the other model parameters are 
optimized. The result is shown in Figure 2. The minimum \ 2 is very small because the degree 
of freedom (hereafter DOF) is and therefore the mass model that reproduces the observables 
perfectly can exist, as is often the case with the analysis of strong lens systems (e.g., Keeton, 
Winn 2003; Pindor et al. 2005). This confirms earlier claims that only image positions and 
fluxes cannot constrain the inner mass slope a very well. 

Figure 3 shows the Ax 2 of the models for each a. Their distribution have no distinct 
feature and therefore local minima of x 2 are expected to exist in a broad range of the parameter 
space. Thus the predicted time delays should also have a broad range of values. In what follows 
we consider only local minima and neglect the distributions of \ 2 around local minima, because 
the differences of time delays between different local minima are typically much larger than the 
uncertainties of time delays around local minima (see also Keeton, Winn 2003). The time 
delays that local minima predict are shown in Figure 4 for a = 0.5, 1.0, 1.5. As shown in Oguri 
et al. (2004, for a = 1.0), all the two of the time delays, for instance the time delay between A 
and B, and C and D, are approximately proportional to each other. We find that this is also 
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Fig. 4. Predictions of the three time delays for different inner slopes: a = 0.5 (red circles), 1.0 (green 
crosses), 1.5 (blue triangles). We show time delay predictions for a group of acceptable models: Each 
point indicates predicted time delays of one of the models. 

the case for a 7^ 1. However, the dispersion is so large — one short time delay does not always 
give the other two short time delays. The distribution in longer time delays is sparse relative 
to that in shorter ones. As increasing a, the distribution shifts to longer time delay space and 
the dispersion becomes larger. There is different trends in the short time delay edges in the 
time delay space for different values of a. Thus, the ratios between time delays are expected 
to vary by the inner mass slope a. 

The tendency of increasing time delays with increasing inner slope a may be ex- 
plained in terms of the famous radial mass index versus the Hubble constant degeneracy 
(Wambsganss, Paczyhski 1994; Oguri et al. 2002; Wucknitz 2002; Kochanek 2002; Oguri, 
Kawano 2003; Kochanek et al. 2005). Wucknitz (2002) showed that modeling a lens galaxy with 
observed time delay as having convergence k oc r^~ 2 results in the degeneracy oc (2-/3). 
More properly, larger f3 leads to smaller h 70 Atij when all the images exist at similar radii. 
The inner mass distribution of the GNFW profile is r~ a , hence its projected mass profile is 
approximately oc r~ a+1 in the case of the GNFW profile. Thus, larger a leads to longer Ai^ 
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Fig. 5. Probability distribution functions of time delay ratios for inner slopes a = 0.5 (red lines). 1.0 
(green lines), 1.5 (blue lines). Their normalizations are arbitrary. 

for the fixed Hubble constant. We confirm that this tendency, at least qualitatively, exists even 
in SDSS J1004+4112 for which the existence of the central galaxy Gl and substructures in 
clusters complicates the total mass distribution. 

To evaluate distributions of predicted time delays, we construct the probability distri- 
bution functions (PDFs) of all the six time delays (At B A, At CB , At CD , Ate a, At ad, At BD ) 
and the ratios between them, by summing up all acceptable models with the weight of 
exp[— (x 2 ~ Xmin)/^]- Again, we neglect the distributions of x 2 around each local minimum, 
since the predicted time delays does not significantly change around each local minimum and 
the uncertainties of time delays are similar for different local minima. Because of the fact that 
width of the valley of the local minimum is not significantly changing, we can exclude the sta- 
tistical weight. Our choice of the maximum x 2 , 11.8, is not problematic because it corresponds 
to more than 99% (A% 2 = 11.3) confidence limit even for DOF = 3. They enable one to show 
that observing the time delays may be able to constrain the mass profile, especially the slope 
a. The PDFs are constructed as 
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fitting parameters are summarized in Table 2. 
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Pij(At\a) oc y^exp 
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S(At - At 



(10) 



where subscript min denotes the minimum \ 2 (f° r each fixed a). In fact, we divide each time 
delay into 100 cells, which all have regular intervals, and calculate the PDF on each cell. The 
PDFs in two and three time delay spaces are defined in the same way. Similarly, the cumulative 
probability distribution function can also be defined as 

P ij {<At\a) = / P i:j (At'\a)d(At'). (11) 

J 

Before going to the PDFs of time delays themselves, we see the PDFs of time delay ratios 
in Figure 5. They slightly broaden but clearly reveal the proportionality seen in Figure 4. The 
PDFs of ratio A^cd/A^ba for a = 1.0 is in good agreement with the value 143 ± 16 of Oguri 
et al. (2004). However, mean values of some ratios are changing significantly with changing a. 
For example, the ratio Atco / At-QA for & — 0.5 and 1.5 are ~ 155 and ~ 125, respectively. The 
ratios A^ad / A^ba and At B D / A^ba depend weakly on the slope — the value of a has relatively 
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Table 2. Fitting Result 





Ci(95%) 


C 2 (95%) 


Ci(50%) 


C 2 (50%) 


Ci(5%) 


C 2 (5%) 


AtBA 


11.8 


11.0 


6.6 


3.2 


0.4 


5.5 


Ai CB 


609.1 


334.3 


330.0 


-1.7 


132.3 


-29.1 


Ai C D 


2232.3 


949.1 


1125.1 


228.6 


173.1 


582.9 


A^CA 


603.7 


359.4 


350.3 


-4.6 


126.6 


-22.9 


At A D 


1512.0 


810.9 


739.1 


280.6 


79.1 


510.3 


At BD 


1532.3 


818.3 


758.6 


260.0 


89.7 


509.1 



much effect between image C and the other images. Thus, from the measurements of the time 
delay between A and B, we can estimate the time delays between A and D, and B and D in 
good accuracy, regardless of the inner slope a. The dependence of image C on the slope a 
corresponds the edges we have seen in Figure 4, to some extent. We note that this figure is 
important to determine the observing strategy of the time delays after measuring the shortest 
time delay A^ba- The figure also implies that determining multiple time delays leads to tighter 
constraints on the inner mass slope a. 

Next we see the PDFs of time delays. Figure 6 shows the cumulative probability of all 
the six time delays. The probability has clear features which we cannot easily see in the time 
delay space (Figure 4), which are summarized as follows: 

1. The models with smaller a have a large fraction of short time delays, and the shorter 
minimum and maximum time delays. For instance, a = 0.5 model has predicted A^ba of ~ 1 — 26 

days, and its cumulative probability reaches 50% at A^ba ~ 8 hjQ days which is nearly 
equal to the minimum value of a = 1.5 model, while a = 1.5 model has predicted A^ba of 
~ 8 — 82 days, and its cumulative probability reaches 50% at A^ba ~ 1 1 hjQ days. 

2. The larger a model has a longer tail of the probability distribution and hence the 
model uncertainties becomes larger with increasing a. 

3. For only A^cb an d A^ca the probability in the short time delay region gives different 
feature — a dependence cannot be almost seen in a region of the short time delays. 

Put another way, increasing a makes the probability distribution to shift to longer time 
delays and to have longer tails. Thus, in principle the measurement of time delays will constrain 
an acceptable range of a. However, only one observed time delay, for example A^ba which is 
the shortest one and likely to be measured most easily, constrains a weakly, given the large 
overlap of the PDFs seen in Figure 6. The exception is the observation of very short or long 
time delay, which gives us the higher or lower limit on the value of a. Because of the third 
feature, it is supposed that constraint from the time delay between image B and C is weak, if 
the time delay is relatively short. 

As described above, the famous mass-slope versus h 70 Atij degeneracy may explain that 
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larger a models give longer time delays. To test this more explicitly, Figure 7 shows time delays 
at which the cumulative PDFs reach 5%, 50%, and 95% for different values of the inner slope 
a. We can fit them with simple linear lines well. Specifically, we fit them as 



We summarize the fitting parameters C\ and C 2 in Table 2. But the time delays do not exactly 
obey the linear law of the scaling relation oc (a — 1). While the discrepancy may reflect the 
complicated mass distribution of the lensing cluster, we still have the qualitative relationship 
between the slope of the mass distribution and time delays. 

4. Constraints on the Radial Slope of the Cluster Mass 

From the discussions in the previous section, it is not obvious how well we can constrain 
the inner slope a from the measurements of time delays. Thus, in this section we present 
a Monte-Carlo analysis for this. First, we assume the value of the inner mass slope a and 
randomly pick up a model from 100 models which are used to construct the PDFs. In other 
word, we assume that one, two, or three of time delays the model predicts are "observed." The 
errors of the "observed" time delays are assumed to be 10% Gaussian because most of observed 
time delays so far have errors of ~ 10%. Then, we calculate the likelihood function 



with Atjj obs and a Ati . = 0.lAt i: , jObs being the value of randomly generated "observed" time delay 
and its error, respectively. Hereafter we consider the three situations: 1) only one time delay 
(A^ba) is measured, 2) two of three time delays (At B A, At CB ) are measured, and 3) all the 
three time delays (At B A, Ai C B, At CD ) are measured. For the cases of multiple measurements 
of time delays, the likelihood becomes multiple-order integral. For each a, we use nine models 
which are randomly chosen from 100 fitted models. We calculate the integral on each cell in 
same way as calculating PDFs. Table 3 shows time delays in units of days, which we 
adopt to constrain the slope a. For instance, when we assume that three time delays of a 
model of (1.25-6), we insert At BAj0bs = 27.4, At CBi0bs = 1353.4, At CD , obs = 4567.2, a AtBA = 2.74, 
CAt CB = 135.34, o"At CD = 456.72 in units of /ifgdays into equation (14), and then from equation 
(13) we calculate the likelihood function for each a with the PDF constructed in §3. 

The results are shown in Figure 8. Each thin and thick solid line show likelihoods 
from each realization and likelihoods averaged over nine realizations. From only one time 
delay measurement, one cannot constrain the inner mass profile a significantly, as expected. 
But by measuring multiple time delays the inner slope can be constrained with good accuracy 



Aty (fij(< Atijla) = 5%, 50%, 95%) = (d + C 2 a)h^dajs. 



(12) 




(13) 
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Fig. 8. The likelihood function £/£ max as a function of the inner slope a. From top to bottom, the 
input a, denoted by short dashed lines, are 0.5, 0.75, 1.0, 1.25, 1.5, respectively. The numbers of time 
delays for constraining the slope are one (A^ba), two (A£ba, A£cb), and three (AtsA, AicB, A£cd) from 
left to right panels. Thin solid lines show the distributions of £/£ max for each realizations (only 5 out 
of 9 realizations are shown, mainly because of the illustrative reason), and thick solid lines indicate the 
distributions after averaging over all 9 realizations. 

(ff Q < 0.25). While each time delay, especially AtcB, gives us little information, the combination 
of time delays leads us to a true value of the slope. We then conclude that longer time delays 
in addition to the shortest time delay are necessary to obtain tight constraints on the inner 
slope a. 

5. Discussion and Conclusions 

We have presented predictions of time delays for the giant quadruple lensed quasar SDSS 
J1004+4112, and investigated the relationship between the time delays and the inner slope of 
the lensing cluster. We adopt a two-component model in which the brightest cluster galaxy 
and the cluster are described by the singular isothermal ellipsoid and the generalized NFW 
profile. We parameterize the inner slope of the cluster component by a such that p oc r~ a in 
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the innermost region. The values of the slope a we have considered are 0.5, 0.75, 1.0, 1.25, and 
1.5, while fixing the scale length to be 40". We have derived a group of mass models that fit 
the observables, and calculated the range of predicted time delays for each a. For observables, 
we have used the data of image positions, fluxes, and the position angle of the galaxy to predict 
the time delays. We have obtained a set of 100 acceptable models for each a, and constructed 
the PDFs of the predicted time delays and the ratios between them. 

We have found that predicted time delays indeed depend on the inner slope, such that 
the steeper inner profiles (larger a) predict longer time delays. All the two of them are approx- 
imately proportional to each other, but the ratios depend slightly on the inner slope, which 
suggests that different a leads to a different structure of the PDFs of multiple time delays. The 
larger a models predict longer minimum and maximum time delays, and have longer tails at 
the higher end of the cumulative probabilities. The time delays At CB and A£ C a give slightly 
different feature in short time delay — the models with different a show almost same distri- 
butions. It is interesting that the model uncertainties resemble the inner slope uncertainty in 
prediction of time delays, and it would be a clue to investigate the model uncertainties. 

To illustrate how well we can constrain the inner slope by adding measurements of time 
delays, we have calculated the likelihood function for a. Figure 8 has shown that the slope 
a is constrained weakly with the measurement of one time delay. However, we also found 
that multiple time delays will results in reasonably strong constraints (a a < 0.25) on the inner 
slope. We note that we assumed observational errors of time delays to be 10%, which may 
be too conservative. For instance, Q 0957+561 has a time delay of ~ 420 days and the error 
is measured to be 1% (Kundic et al. 1997; Oscoz et al. 2001). If this level of errors can be 
achieved for SDSS J1004+4112, we will be able to determine the inner slope more tightly. 

Our model predictions offer useful guidance for photometric monitoring of this lens 
system to determine the time delays. We found that 95% of models of a = 0.5 — 1.5 have 
the predicted time delays of A^ba ^5 28, AtcB ^5 1400, and A^cd ^ 3700 in units of h^Q days. 
The first time delays is similar to those in any other lensed quasars, and therefore it can be 
measured easily. Indeed, the preliminary detection of the time delay between B and A has 
been made to be ~ 25 days (C. S. Kochanek, private communication). The second shortest 
time delay A^cb is ~ 30 — 40 times larger than AteA; thus assuming the measured AteA we 
predict A^cb = 750 — 1000 days. This is somewhat longer than observed time delays in any 
other systems, but is not impossible to measure. The longest time delays need monitoring more 
than ~ 10 years, making the measurement quite challenging. 

The mass profile of the lensing cluster can be constrained much better if we combine 
the time delay measurements with other observations. For instance we may add constraints 
from multiply imaged galaxies behind the cluster (Sharon et al. 2005). The multiple arcs 
that have different redshifts (and therefore different Einstein ring radii) are an independent 
tool to constrain the cluster mass profile, especially the radial mass slope a. Other such 
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observations includes X-ray measurements of intra-cluster medium and the measurement of the 
velocity dispersion of the brightest cluster galaxy (eq. [8]). By combining these complementary 
information, we will be able to reveal the detailed distribution of the dark matter of the lensing 
cluster in a robust manner. 

We are grateful to Takahiko Matsubara, Naohisa Inada, and Chris Kochanek for discus- 
sions. 

References 

Bourassa, R. R., & Kantowski, R. 1975, ApJ, 195, 13 

Bourassa, R. R., & Kantowski, R. 1976, ApJ, 205, 674 

Bourassa, R. R., Kantowski, R., & Norton, T. D. 1973, ApJ, 185, 747 

Broadhurst, T., et al. 2005, ApJ, 621, 53 

Chiba, M. 2002, ApJ, 565, 17 

de Blok, W. J. G., & Bosnia, A. 2002, A&A, 385, 816 
Diemand, J., Moore, B., & Stadel, J. 2004, MNRAS, 353, 624 
Freedman, W. L., et al. 2001, ApJ, 553, 47 
Fukushige, T., & Makino, J. 1997, ApJ, 477, L9 
Fukushige, T., & Makino, J. 2001, ApJ, 557, 533 
Fukushige, T., & Makino, J. 2003, ApJ, 588, 674 
Fukushige, T., Kawai, A., & Makino, J. 2004, ApJ, 606, 625 

Gentile, G., Salucci, P., Klein, U., Vergani, D., & Kalberla, P. 2004, MNRAS, 351, 903 
Ghigna, S., Moore, B., Governato, F., Lake, G., Quinn, T., & Stadel, J. 2000, ApJ, 544, 616 
Hayashi, E., et al. 2004, MNRAS, 355, 794 

Hennawi, J. F., Dalai, N., Bode, P., & Ostriker, J. P. 2005, ApJ, submitted (astro-ph/0506171) 

Inada, N., et al. 2003, Nature, 426, 810 

Inada, N., et al. 2005, PAS J, 57, L7 

Jing, Y. P., & Suto, Y. 2000, ApJ, 529, L69 

Kaiser, N., & Squire, G. 1993, ApJ, 404, 441 

Kawano, Y., Oguri, M., Matsubara, T., & Ikeuchi, S. 2004, PASJ, 56, 253 

Keeton, C. R. 2001, preprint (astro-ph/0102340) 

Keeton, C. R., Kochanek, C. S., & Falco, E. E. 1998, ApJ, 509, 561 

Keeton, C. R., Kochanek, C. S., & Seljak, U. 1997, ApJ, 482, 604 

Keeton, C. R., & Winn, J. N. 2003, ApJ, 590, 39 

Kochanek, C. S. 2002, ApJ, 578, 25 

Kochanek, C. S., & Dalai, N. 2004, ApJ, 610, 69 

Kochanek, C. S., Morgan, N. D., Falco, E. E., McLeod, B. A., Winn, J. N., Dembicky, J., & Ketzeback, 

B. 2005, ApJ, submitted (astro-ph/0508070) 
Kundic, T., et al. 1997, ApJ, 482, 75 
Mao, S., & Schneider, P. 1998, MNRAS, 295, 587 



18 



Moore, B., Quinn, T., Governato, F., Stadel, J., & Lake, G. 1999, MNRAS, 310, 1147 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 462, 563 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

Oguri, M., et al. 2004, ApJ, 605, 78 

Oguri, M., & Kawano, Y. 2003, MNRAS, 338, L25 

Oguri, M., & Keeton, C. R. 2004, ApJ, 610, 663 

Oguri, M., Taruya, A., Suto, Y., & Turner, E. L. 2002, ApJ, 568, 488 

Oscoz, A., et al. 2001, ApJ, 552, 81 

Pindor, B., et al. 2005, AJ, in press (astro-ph/0509296) 

Power, C, Navarro, J. F., Jenkins, A., Frenk, C. S., White, S. D. M., Springel, V., Stadel, J., & 

Quinn, T. 2003, MNRAS, 338, 14 
Reed, D., Governato, F., Verde, L., Gardner, J., Quinn, T., Stadel, J., Merritt, D., &: Lake, G. 2005, 

MNRAS, 357, 82 
Refsdal, S. 1964, MNRAS, 128, 307 
Saha, P., & Williams, L. L. R. 1997, MNRAS, 292, 148 
Salucci, P., & Burkert, A. 2000, ApJ, 537, L9 

Sato, S., Akimoto, F., Furuzawa, A., Tawara, Y., Watanabe, M., & Kumai, Y. 2000, ApJ, 537, L73 
Schneider, P. 2005, astro-ph/0509252 

Schneider, P., Ehlers, J., &: Falco, E. E. 1992, Gravitational Lenses (New York: Springer) 
Sharon, K., et al. 2005, ApJ, 629, L73 

Spergel, D. N., & Steinhardt P. J. 2000, Phys. Rev. Lett., 84, 3760 
Spergel, D. N., et al. 2003, ApJS, 148, 175 

Swaters, R. A., Madore, B. F., van den Bosch, F. C, & Balcells, M. 2003, ApJ, 583, 732 

Tasitsiomi, A., Kravtsov, A. V., Gottlober, S., & Klypin, A. A. 2004, ApJ, 607, 125 

Walsh, D., Carswell, R. F., & Weymann, R. J. 1979, Nature, 279, 381 

Wambsganss, J., & Paczyhski, B. 1994, AJ, 108, 1156 

Weinberg, M. D., & Katz, N. 2002, ApJ, 580, 627 

Weymann, R. J., et al. 1980, Nature, 285, 641 

Williams, L. L. R., & Saha, P. 2004, 128, 2631 

Wucknitz, O. 2002, MNRAS, 332, 951 

York, D. G., et al. 2000, AJ, 120, 1579 



19 



Table 3. Input Time Delays in Likelihood Analysis 
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Table 3. (Continued.) 
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* Name of data arc expressed as (input a-data number), 
t Time delays in units of h^ . The observational errors 
are assumed to be 10%. 
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